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Abstract 

This paper studies tlie quantum dynamics of a charged particle in a 2D square lattice, under 
the influence of electric and magnetic flelds, the former being aligned with one of the lattice axes 
and the latter perpendicular to the lattice plane. While in free space these dynamics consist of 
uniform motions in the direction orthogonal to the electric fleld vector, we flnd that, in a lattice, 
this directed drift takes place only for specific initial conditions and for electric field magnitudes 
smaller than a critical value. Otherwise, the quantum wave-packet spreads ballistically in both 
directions orthogonal to the electric field. We quantify this ballistic spreading and identify the 
subspace of initial conditions insuring directed transport with the drift velocity. We also describe 
the effect of disorder in the system. 
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I. INTRODUCTION 



Quantum transport in periodic potentials is a topic of permanent interest since the early 
days of quantum mechanics. Recently this topic has attracted a renewed interest thanks to 
the experiments with cold atoms in optical lattices, where the quantum dynamics can be 
observed both in real and in momentum space. In particular, during the last decade much 
attention has been devoted to the phenomenon of Bloch oscillations of cold atoms subject to 
a static (for example, gravitational) field (see e.g. Ref. [l-4| to cite just a few of the relevant 
papers). These new studies have also shed additional light on the old problem of electron 
transport in a solid crystal induced by an external electric field 5|-l3]. 

Present research in cold atoms physics is also focused on the problem of generating 
synthetic magnetic fields, which could impart a Lorenz-like force to otherwise neutral atoms 



in motion 



This opens an interesting perspective for studying the quantum Hall effect 



with cold atoms in 2D optical lattices 



12| . A preliminary but necessary step in this direction 



is the analysis of the Bloch dynamics of cold atoms in the presence of an artificial magnetic 



field 13]. In this paper we study these dynamics in the tight-binding approximation, where 
the single-particle Hamiltonian of an atom in a 2D lattice has a relatively simple form. In 
spite of this simple form, we find a surprisingly rich variety of behaviors, which range from 
the ballistic spreading of the wave-packet to directed transport. 

The structure of the paper is as follows. In Sec. [TTl we introduce the model and perform 
a preliminary analysis, following the standard route which leads to the Harper equation 



161] and, hence, is referred to as the Harper approach. The Harper approach reduces the 



original 2D problem to a ID problem, which on the one hand is simpler to describe, but on 
the other hand it limits the class of initial conditions that can be considered. Therefore, in 
Sec. mil we tackle the problem using the more powerful formalism of Landau-Stark states, 
which are the eigenstates of a charged particle in a lattice subject to both magnetic and 
electric fields. In some sense this approach is opposite to the Harper approach: here one 
begins with a particle in an electric field (i.e., with the problem of Bloch oscillations) and 
then introduces a magnetic field. The Landau-Stark states approach is fully two dimensional 
and, hence, imposes no limitations on the class of initial conditions that can be considered. 
In addition, it allows us to construct 2D localized wave packets, which propagate across the 
lattice without changing their shapes. 
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Finally, in Sec. IIVI we discuss the transport in the presence of on-site disorder. Here 
we consider generic initial conditions in the form of a wide incoherent wave packet. In this 
case, for vanishing disorder, the wave-packet spreads ballistically in the direction orthogonal 
to the electric field. We study the modifications of this ballistic spreading caused by the 
presence of on-site disorder. In the case of null magnetic field this problem was considered in 



the recent paper [17|. Our analysis unveils the differences between the ID Harper approach 
of Sec. [IT] and the full 2D approach of Sec. Illlt that, to the contrary, predict the same rate of 
ballistic spreading when disorder is absent or negligible. This provides us with a significant 
new model for the problem of Anderson localization, where the effective dimensionality is 
larger than one but smaller than two. 

II. THE HAMILTONIAN MODEL AND ITS APPROXIMATIONS 

In this section we introduce the physical model under investigation and we discuss two 
common approximations. We also derive a first dynamical phenomenon: the presence of a 
drift of the particle/wave-packet with constant velocity. Let us therefore consider a quantum 
particle of mass M and charge e in a 2D square lattice of side d in the x-y plane, created 
by a periodic potential V{x, y). The particle is also under the action of an in-plane electric 
field, aligned with the y axis and of a magnetic field normal to the x-y plane. This particle 
is described by the Hamiltonian 

H = ^(y^--^P^ +V{x,y) , V{x + ld,y + md) = V{x,y) , (1) 

where A is the vector potential. A convenient choice of the vector potential, corresponding 
to the field configuration described above is 

A = fi(0,x,0) + F(0,ct,0) , (2) 

where B and F are the magnetic and electric field magnitudes, respectively, and where 
vectors are described in line notation. It is easy to check that the Hamiltonian ([1]) commutes 
with the generalized translation operator = exp(— z27ra?//(i)T^., where is the usual 
translation operator and where 

eBd' 

The parameter a defines the so-called magnetic period X = d/a, which has the dimensions of 
a length and will play an important role in the following. A rigorous analysis of the motion 



generated by the Hamiltonian ([T]) is a difficult, unsolved problem. Therefore, one deals with 
this problem by introducing suitable approximations. 



A. The effective mass approximation 

The simplest approximation of the Hamiltonian ([1]) starts by defining the effective mass 
M* as the curvature of the ground Bloch band of ([1]) at A = 0: 

M* = {M^Myf'\ M.,. = ^^, (4) 

where E = E{kx, Ky) is the dispersion relation. The approximate Hamiltonian then reads 

For vanishing electric field the eigenfunctions of (|5]) are Landau oscillators with energy 
spectrum En = huc{n + 1/2). Here Uc = eB/cM* is the cyclotron frequency. It appears 
from eq. ([5]) that a non-zero value of F shifts the origin of a Landau oscillator linearly in 
time. Thus, the particle is transported in the direction orthogonal to the electric field with 
a drift velocity v* = cF/B. 

Although the effective mass approximation correctly predicts the drift velocity, it com- 
pletely ignores the lattice discreteness, which is important in many aspects. To take the 
lattice discreteness into account one usually introduces a tight-binding approximation to 
the original Hamiltonian. 

B. The tight-binding approximation 

The tight-binding approximation amounts to the following ansatz for the wave function 
of the system, 

that is written in terms of two-dimensional Wannier states Wi^m{x,y) = \l,m). In this 
approximation one also neglects coupling other than between nearest neighbors pairs of 
states, so to approximate the Hamiltonian (IT]|3]) in the form 

= -y E (1^ + 1' ^) ^1 + h.c.) - ^ E (l^. + 1) (m|e^(2™;-a.,i) ^ ^^^^^ ^ 

l,m l,m 
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where ub = edF/h is the Bloch frequency. This latter is the second characteristic frequency 
of the system, which is overlooked by the effective mass approximation. The hopping matrix 
elements J^^ and Jy in the Hamiltonian ([7]) can be related via eq. to the parameters of 
the original system as J^^y = h'^ /(PM* y. In what follows, to simplify equations, we shall 
set e, c, h, M and d to unity. The magnetic period is then given by A = the Bloch 
frequency becomes ub = F, the cyclotron frequency 

LJc = 2na{J,JyY/^ , (8) 

and the drift velocity 

V* = F/27ia . (9) 

The Hamiltonian (jTj) with four parameters Jx, Jy, a, and ub = F defines the model being 
analyzed in the following sections. More precisely, we are interested in wave-packet dynamics 
of the system ([7]), which is governed by the Schrodinger equation for the amplitudes ipi^mit)'- 

Remark that the discrete index / refers to the x direction, while m labels the y direction, 
parallel to the electric field. 



C. Semiclassical analysis of a one-dimensional reduced Hamiltonian 

We begin our analysis of the system dynamics with considering a class of initial conditions 
with uniform probability density along the y direction. This leads to a one-dimensional 
reduction of the quantum problem. In fact, in this case we may use the substitution 

gi/tm 271" 
^Pl^^it) = -^bi{t) , K = —k, (11) 



which reduces the Schrodinger equation (11 01) to the following ID equation for the amplitudes 
hiit): 

ibi = -y (&i+i + k-i) - Jy cos{27ral + k- Ft)bi . (12) 

The substitution ( ITTj) assumes periodic boundary conditions, of period Ly, that eventually 
tends to infinity. Equation f|T2|) can be rewritten in the more familiar form of the Schrodinger 
equation for the driven Harper Hamiltonian 

Hid = — Jx cosp — Jy cos(27rQ;a; + k — Ft) . (13) 
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FIG. 1: (color online). Stroboscopic map of the driven Harper for = Jy = 1 and F = 0, a small 
F <^ Fcr and moderate F below and above Fcr- 



In this paper we shall mostly limit ourselves to the case a <^ 1, which can be considered as a 
semiclassical limit for Harper-like Hamiltonians 18|]. Thus, we can appeal to the quantum- 
classical correspondence to facilitate the physical understanding of the problem. Letting the 
operators x' = 2iTax and p' = —i2Tia d/dx' be considered as classical canonical variables 
yields the the classical driven Harper Hamiltonian 



H, 



d 



-Jx cos(p') — Jy cos(x' — Ft) . 



(14) 



Figure [T] shows the Poincare' surfaces of section for the classical motions. Phase trajec- 
tories near the stationary point {p',x') = (0,0) in the panel (a), which refers to the case 
F = 0, encircle this point with the cyclotron frequency (|8]) and are thus associated with 
the Landau states in the effective mass approximation. For weak field intensity F, these 
trajectories are captured into the nonlinear resonance, which is seen as the stability island 
in Fig. (H^b-c), and are transported with the drift velocity v*. The stability island, however, 
shrinks with increasing F and completely disappears when > Fcr, where 

Fcr = 27taJ^ . (15) 
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Thus the drift is possible only if |F| < Fcr- Translated into the quantum language this 
means that the Landau state can be transported only if the condition f lTSj) is satisfied. To 
support this statement, in Fig. |2]we show the dynamics of a localized wave packet associated 
with the ground Landau state in the effective mass approximation, given by the initial 
condition = 0) ~ exp(— 7ra/^). A nice, directed transport is observed only in a weak 
field regime, while in the strong field regime the wave packet spreads almost symmetrically 
both directions. Remark that these results refer to the one-dimensional model f|T2|) . We will 
describe wave-packet spreading in more detail in Sec. HVAl 

We also note that eq. ffT^ removes a seemingly illogical consequence of eq. that 
predicts an infinite drift velocity for null magnetic field. As a matter of facts, when the 
magnetic field intensity B tends to zero, the transporting island in Figure [T] disappears. 
Thus, eq. iQ for the drift current implicitly assumes that F < 2T[aJx- 

To conclude this section we note that directed transport in the weak field regime (i.e., 
the drift) takes place not only for states associated with trajectories near = (0,0), 

but also for states associated with trajectories near {p',x') = (vr,7r) in Fig. [l](a). These 
are high-energy counterparts of the low-energy (Landau) states of the Harper Hamiltonian, 
related by the transformation 5^*^^^ = [—iybf°^\ 



III. QUANTUM ANALYSIS OF THE 2D SYSTEM 

We now move to the analysis of the two-dimensional tight-binding Hamiltonian. By using 
the Kramers-Hennenberger transformation, ipi^mit) '^i,m{t) e^p{—iFmt), we can reduce 
the Hamiltonian (jT]) to the form 

This can also be seen as the tight-binding approximation to the Hamiltonian 

^=^(p-^a) +V{x,y) + eFy, 

where the vector potential A = B{0,x,0) is used instead of that in eq. ([2]). They describe 
the same physical situation. 
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FIG. 2: Space-time plot of the dynamics of a wave packet initially localized about the origin for 
Jx = Jy = ^, n = 0, a = 1/20 and F = 0,0.5,0.1 (left, center, right). To increase visibility, the 
packets for F = and F = 0.1 have been shifted by ^20 lattice sites. Wave-functions amplitudes 
|6/(f)P are reported in grey scale, black maximum. The time axis is measured in units of the Bloch 
period Tg = 2ir/uJB and we have arbitrarily set Tg = 2it for null electric field. For i*" = 0.1 the 
particle is captured into the nonlinear resonance and transported across the lattice, traveling one 
magnetic period A in one Bloch period Tg. 

A. Landau-Stark states 

One gets important insight into the dynamics by analyzing the energy spectrum of the 
Hamiltonian f|T6l) . (Htiy-ip)i^rn = Eipi^^. By assuming periodicity in the x direction, of an 
arbitrary (large) period L^., and eventually letting tend to infinity, we can write the 
eigenf unctions in the form 



where k is free to vary in the interval [0, 27r]. [Please notice the difference between eqs. ( IT711 
and ( ITT]) .] This reduces the spectral problem to the one- dimensional problem: 




,m 




(17) 



(cm+i + Cm-i) + [Fm — J.J. cos(27ram — n)\cm = Ec, 



(18) 
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FIG. 3: (color online). A fragment of the energy spectrum of (|16p for J^. = Jj, = 1, a = 1/10, and 
F =1 (left) and F = 0.3 (right). 

Labeling the eigen-solutions of f|T8|) by the discrete index z/, so that E = E^{k,) and Cm = 
c^{i^) and inserting these values in eq. (fT7|) . we construct the eigenfunctions of (fT6|) in the 
form of the Landau-Stark states: \^u,t^ = Y,i,m'4'trn'\h^) ■ 

For = the spectrum of (ITSl) consists of a ladder of Bloch bands of zero width, 
separated by the Stark energy F. For ^ these bands acquire a finite width ~ 2J2, 
and for any fixed k, the spectrum is a modulated Wannier-Stark ladder. A fragment of the 
spectrum covering three modulation periods is shown in Fig. [3] for two values of the field 
intensity F. Note that for small F a pattern of straight lines emerges in the picture. The 
eigenstates belonging to this lines (which always appear in pairs) are associated with two 
transporting islands in the classical approach. We shall discuss them in more details in the 
next subsection. 

It follows from eq. (fT7|) that Landau-Stark states are extended states in the x direction. 
Moreover, they are current states, i.e. the mean value of the current operator 0, 



is non-zero for almost all Landau-Stark states [^v^k)- This value, Vy{K), can be calculated 




(19) 
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FIG. 4: (color online). The mean current (I19p carried by Landau-Stark states \^i/,k) with k = 0.1. 
Parameters are Jx = Jy = 1, F = 0.3, a = 1/10 (upper panel) and a = 1/10.1417 (lower panel). 
The dashed line corresponds to the drift velocity v* . 

via eq. f|T71) : 

v^{k) = (^^,«|{)|^^,«) = |c:;(K)|^sin(27ram - k) . (20) 

m 

The upper panel in Fig. H] shows the mean current for a set of one-hundred states \'^u,k.) 
with K, = 0.1 for the same parameters of Fig. [3](b). It is seen that the Landau-Stark 
states carry currents of different magnitude and sign. However, completeness of the basis of 
eigenstates implies the sum rule 

J2vA^) = 0. (21) 

V 

A remark about commensurably of the magnetic and lattice periods is in turn. In Fig. [3] 
the spectrum Eyi^n) is periodic along the energy axis because a is rational. This periodicity is 
replaced by quasi-periodicity when a is irrational. As a consequence, the same happens also 
for fi/(/t), as can be seen clearly in Fig.Hl^b). However, Landau-Stark states are still extended 
functions in the x direction and the summation rule fl2T]) also holds. Thus, irrationality of 
the value a does not change properties of the system in a crucial way. This differs from the 
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case of vanishing electric field, where the eigenstates are localized functions for irrational a 

if J. > jy m. 



B. Transporting states 

In this subsection we describe a procedure for constructing localized wave packets, which 
propagate across the lattice with the drift velocity velocity v* in eq. ([9]) without changing 
their shapes. 

First of all we note that for \F\ < F^r the potential energy term in eq. f llSp . V{m) = 
Fm — JxCos{27Tam — k), has local minima and maxima. These local extrema can support 
a number of well localized states. The exact number of these states depends on the system 
parameters and coincides with number of straight lines in the energy spectrum. In first 
order approximation they can be found by expanding V{m) in the vicinity of the extremum. 
Without loss of generality, let us consider the minimum at m = 0. This expansion reduces 
(ITS]) to the following equation. 



(the constant term is omitted), which is a Mathieu-type equation whose properties are 
well-known. As an example. Fig. |5] compares the localized state, which is calculated using 
the exact equation (ITSll with that following from the approximate equation ( l22l) for the 
parameters of Fig. [3]^b). The two states well coincide with the exception of ten points in the 
K axis, in coincidence of which [see Fig. [3t^b)] the spectrum shows avoided crossings. 

Ignoring avoided crossings [i.e., using eq. ( 122|) instead of eq. ( !T8|) ] and tracking the energy 
Ey{K) along the straight lines in Fig. El^b) we can construct a localized state by integrating 
extended Landau-Stark states over the quasimomentum. 



It follows from eq. (122|) that the time evolution of this state is defined by the trivial shift. 



(cm+i + Cm_i) + Fm H — -{2Tiam — KY Cm = Fa, 



■m 



(22) 




(23) 



'4'l,m{t) = iJl-v't,m{0) , 



(24) 



where v* is the drift velocity. Needless to say, transporting states can be constructed for any 
local minimum or maximum of ^(m), i.e. for any straight line in the spectrum. Examples of 
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FIG. 5: Absolute values of the expansion coefficients Cm, calculated on the basis of the exact 
equation (fT8|) and the approximate equation (p2]) . Parameters are Jx = Jy = 1, F = 0.3, and 
a = 1/10. 
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FIG. 6: Grey scale representation of transporting states. Parameters are Jx = Jy = 1, F = 0.1, 
and a = 1/10 (left) and a = 1/20 (right). 

transporting states for F = 0.1 and two different values of a are given Fig. El It is seen that 
the states ( 123|) are well localized in both spatial directions, the characteristic localization 
length along the y axis being defined by the magnetic period X = 1/a. 

We mentioned above that transporting states are associated with transporting islands 
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FIG. 7: (color online). A fragment of the energy spectrum of the system (left) and the 'exotic' 
transporting state (right), which transports particles in the counter- intuitive direction. The system 
parameters are Jx = Jy = 1, F = 0.1 and a = 1/(2.2). 

in the semiclassical approach, which is justified if \a\ <^ 1. An advantage of the Landau- 
Stark states formalism is that one can make predictions also for |a| ~ 1/2, i.e., in the deep 
quantum regime. (Without a loss of generality one may consider —1/2 < a < 1/2.) In 
particular, for a ~ 1/2 we have found exotic transporting states, which have no classical 
analogue. An example is given in the left panel of Fig. [TJ which shows the band spectrum 
Ejj{k) for a = 1/(2.2). Using the procedure described above one can construct localized 
transporting states for any straight line in the spectrum, see right panel in Fig. [71 Note 
that, while in Fig. Et^b) the slope of the straight lines is positive, in Fig. El^a) it is negative: 
therefore the depicted localized state transports quantum particles in the opposite direction 
to what (classically) expected by the electric-magnetic field geometry. 

In closing this section we need to stress again that we ignored avoided crossings when 
constructing the transporting states. The presence of these avoided crossings leads to back- 
scattering effects: the traveling packet in eq. flMj) emits tiny packets, which propagate in the 
opposite direction. Thus, strictly speaking, the perfect transport of a localized wave packet 
in a lattice is impossible. 
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IV. WAVE PACKET SPREADING 



The band spectrum of the Landau-Stark states and the summation rule (12T]) imply that 
the time evolution of a generic localized wave packet is a ballistic spreading along both the 
positive and the negative x direction, with the obvious exception of specific initial conditions 
such as the transporting states. In this section we study wave-packet spreading in some 
detail, extending the analysis also to the case of disordered lattices. To construct a generic 
initial wave packet, we shall superimpose random phases to a wide Gaussian envelope of 
width 1/a. Note that by averaging the results over such random phases we can mimic 

the dynamics of an incoherent wave packet. 

A. One-dimensional system 

To have a reference model for the phenomena under investigation, we first analyze the 
dynamics in the one-dimensional system (fT3|) . The simplest situation is obtained for large 
F ^ FcTi where the classical phase-space trajectories are nearly straight lines parallel to 
the X axis. In this case we can safely neglect the potential energy term and hence assume 
that the spreading motion coincides with that in the absence of driving. For an incoherent 
wave packet this implies that the square root of the second spatial moment, cr(t), is 



Equation (125|) gives an upper boundary to the rate of ballistic spreading of a localized wave 
packet. Decreasing F the rate of spreading also diminishes. At the same time, the spreading 
becomes asymmetric, a fact especially evident for F < Fcr- In this regime a large part of the 
packet moves in the positive direction with the drift velocity v*, while a smaller part moves 
in the negative direction with velocity |f | > v*. Note that this asymmetric spreading does 
not imply a preferable transport in either direction, since the first moment of the position 
of the packet is null. Fig. [10] below provides a numerical confirmation of this fact. It is 
obtained by averaging over the random phases of the initial wave packet. 

Next we study the effect of on-site disorder, which we mimic by adding to the ID Hamil- 
tonian in eq. ( IT3l) random on-site energies e„, |e„| < e/2. The dashed lines in Fig. [8] show 
a{t) versus t for increasing values of the disorder, at fixed F = 3. Here, and in the following, 
we also perform an average over disorder realizations. As expected, the ballistic spreading 




(25) 
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FIG. 8: (color online). Wave-packet width as a function of time for = Jy = 1, a = 1/10, 
F = 3 and increasing strength of the disorder e = 0,0.1,0.2,03,0.4,0.5, 1, (color labelled as blue, 
red, magenta, green, blue, red, magenta, respectively) from top to bottom. The continuous lines 
correspond to the 2D model, while the dashed lines to the one-dimensional model. 

is suppressed by the disorder. Moreover, for large times the spreading is completely stopped 
by the Anderson localization. In this case the spatial distribution of the average packet, 
P{l,t) = (iV'i(^)P) approaches the exponential distribution. 



which is a hallmark of Anderson localization. Let us also note that the value F = 3 chosen 
in Fig. [8] is large enough to approach the free expansion regime, eq. fl25|) . for e = 0. Thus 
we expect that for e 7^ the functional dependence of the localization length £ = £(e) in 
eq. (l26l) might coincide with that of the standard Anderson model. 

The case F < Fcr, shown in Fig. [9] (dashed lines), appears to be more complicated. Here 
the regime of ballistic spreading changes to the regime of Anderson localization through 
an intermediate regime of diffusive spreading, where the spatial distribution function of the 
wave packet grows approximately like the diffusion law, P{l,t) ~ exp{—P/Dt). At F = 0.3 
this regime is clearly observed in the numerical simulations for disorder strengths around 



P(/,t)~exp(-|/|/£), 



(26) 
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FIG. 9: (color online). The same as in the previous figure, now for F = 0.3. 



e ~ 0.5. To see that this is actually an intermediate regime, [i.e. to see a sign of saturation) 
we will perform a more sophisticated analysis, in the next section. 

B. Two-dimensional system 

As in the ID case we consider a Gaussian incoherent packet as initial condition for the 
motion, of width ctq and ctq in the two directions; the former is chosen to coincide with 
the width of the initial condition in the one-dimensional model. For vanishing disorder the 
packet spreads ballistically in the x direction but remains localized in the y direction. The 
Landau-Stark localization length Cy in the y direction is in fact: 



The estimate ( 1271) is obtained on the basis of the ID Wannier-Stark problem (ITSl) and 
assumes a <^ 1 and Jx = Jy (note however [20:]). It follows from eq. (1271) that the localization 
length increases from unity to one magnetic period when F is decreased. We also mention 
that in the under-critical regime, F < F^r, the estimate Cy ~ 1/a gives the maximal 
localization length, while the minimal localization length scales as Cy ~ 1/y/a. Thus in 




max(l,2Jj^/F) , F > F, 
1/a , F < F, 



cr 



cr 



(27) 
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FIG. 10: (color online). Upper panel: Population of the lattice sites at the end of numerical 
simulation. Lower panel compares integrated population of the lattice sites along the x direction 
(solid line) with prediction of the ID model (fTO]l . dashed line. 

this regime different Landau-States have different localization lengths £y, in the interval 

l/y/a < Cy < 1/ a. 

We can now show the results of the analysis of the full, two-dimensional system (fT6|l . 
The upper panel in Fig. [TO] depicts the averaged 2D wave-packet (|'?/'/,m(^)P) at t = 2000 for 
the parameters of Fig. [3](b) (note the different scale of the y and x axis). It is seen that the 
packets splits in several sub-packets, of which the rightmost moves with the drift velocity. 
This sub-packet is associated with the straight lines in the energy spectrum, i.e., with the 
transporting states. At the opposite end, the leftmost sub-packet is associated with the part 
of the spectrum where the energy bands have minimal negative derivative. Intermediate 
packets are also observed, so that one can easily predict speed and position of wave-packets 
by analyzing the energy spectrum. 

The lower panel in Fig. [10] compares the projected wave-packet in the x direction, 
Hm I'^z.mP, with the corresponding distribution predicted by the ID model ffT^ . A rea- 
sonable correspondence is noticed, which strengthen the physical significance of the one- 
dimensional reduction. This correspondence is even more precise if one compares integrated 
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characteristics like the wave-packet widths, which we define as 

^W = (E^'|cKi)r)^ and a{t)= (Y^Pl^P, 



l,m I 

I Lm 



in the ID and 2D cases, respectively. This comparison is done in Fig. [H] and Fig. [Hlfor F = 3 
and F = 0.3, respectively. It is seen that the values for e = practically coincide. To the 
contrary, for non-zero amplitude of the disorder, the data of the one-dimensional model are 
consistently smaller than the other, supporting the almost obvious observation that in the 
2D model Anderson localization is "less effective". To turn this observation into rigorous 
theory is the goal of further investigation. 

Here, we only perform a simple analysis. From the data of Fig. [H]and Fig. [Hlone might 
think that anomalous diffusion is taking place in the dynamics of the wave-packet, as in 
well known examples of quantum intermittency (see [2j] and references therein). Indeed, 
this is not the case. In fact, if we compute the local exponent z/(t) as the slope of the curve 
log(cT^(t)) versus logt: 

''^>~ dlogt ' 

we obtain Fig. [TTl that sums up the data for F = 0.3 and F = 3.0. While slopes for e = 
tend clearly to the value = 2, indicating linear motion, the other curves, after reaching a 
maximum, show a clear tendency to decrease. (The initial increase can be easily explained 
as an effect of the wide amplitude of the initial wave-packet.) The successive decrease of the 
local exponent z/(t) is the signature of Anderson localization. In turn, this can be justified 
by noting that Stark localization in the y direction, eq. f l27|) . renders the system quasi one- 
dimensional, effectively bounding the motion to a strip, the larger the value of the electric 
field amplitude F, the narrower this strip. 



V. CONCLUSIONS 



We have studied the wave-packet dynamics of a quantum particle in a square lattice in the 
presence of (real or artificial) electric and magnetic fields, in the tight-binding approximation. 
In this approximation the Hamiltonian of this system contains five physical parameters: the 
hopping matrix elements J^^y, the components of the electric field vector F^^y, and the 
magnetic field fiux quanta through a unit cell a, (with |a| < 1/2). In this paper we have 
considered mainly the case we do not expect significant differences to 
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FIG. 11: (color online). Local exponents v{t) in the growth of the second moment of wave-packets 
as a function of time, computed on the data of Fig. [8] and Fig. [9l Labels: F = 0.3, green hnes, 
asterisks; F = 2>. red hne, crosses. In both cases e = 0, .1, .2, .3, .4 (top to bottom curves). The 
horizontal lines at u = 1 and u = 2 mark the linear and quadratic increase of <t^. 

appear in the more general situation. More importantly, we have assumed that the electric 
field is parallel to one of the lattice axes which, in fact, is a rather specific situation. In 
spite of this limitation the system is found to display very rich dynamics, varying from 
ballistic spreading to directed transport. A laboratory application of this model can consist 
of cold atoms in a 2D optical lattice, subject to artificial electric and magnetic fields. In 
this experiment, typical initial conditions correspond to a coherent (a BEC of atoms) or 
incoherent (a thermal cloud of atoms) 2D wave packet, localized over several lattice sites. 
Our findings apply directly to this experimental situation. 

We have discovered two main regimes for wave-packet dynamics under these geometric 
conditions. These regimes are determined by the relation between two dimensionless pa- 
rameters - the parameter a, which is proportional to magnetic field magnitude and the 
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parameter /3, which is proportional to the electric field magnitude and is defined as the ratio 
of the Stark energy to the hopping matrix element. 

For vanishing magnetic field, a = 0, the electric field induces Bloch oscillations of the 
wave packet in the field direction, while in the orthogonal direction the wave packet spreads 
ballistically with a rate defined by J. Under the condition of Stark localization, /3 ^ 1, this 
regime is also valid if a 7^ 0. Note that the condition /3 ^ 1 implies that the amplitude of 
Bloch oscillations is smaller than one lattice site. 

In the weak field regime, /3 <^ 1, the dynamics are richer and depend on the relative 
size between the parameters a and /3. Namely, if /3 > 27ra the wave-packet still spreads 
ballistically, although the spreading becomes asymmetric. A qualitative change happens if 
P < l-na. Note that in dimensional units this condition reads < Wc, where iob and lOc are 
the Bloch and cyclotron frequencies, respectively. Here, the wave packets splits into several 
sub-packets, of which the right-most moves with the drift velocity v* . This sub-packet is 
supported by the transporting states, the explicit form of which is given in Sec. 1111 B[ 

We have explained the observed dynamical regimes by analyzing the energy spectrum of 
the system, which can be easily calculated numerically. This spectrum has a band structure 
with specific band patterns, sensitive to the system parameters. Analyzing these patterns one 
can make predictions about the wave-packet dynamics even before performing the demanding 
simulation of the 2D Schrodinger equation. 

In this work we have also derived and tested a ID approximation to the original 2D 
Schrodinger equation, which leads to the driven Harper Hamiltonian ( |T3l) . This latter is 
interesting in its own, because the driven Harper system is more easly realized in the lab 
than 2D lattices with artificial electric and magnetic fields 2l|. We have shown that the ID 
model is capable to capture some important features of the 2D wave-packet dynamics. 

Finally we have studied the effect of on-site disorder on the dynamics 22| . It is found that 
disorder changes both the ballistic and transporting regimes into the regime of Anderson 
localization, which can be preceded by a diffusive regime. The observed diffusive spreading 
of the wave packet together with asymptotic Anderson localization indicates that on-site 
disorder converts the extended Landau-Stark states (which are the system eigenstates for 
vanishing disorder) into localized states with a non-trivial scaling law. This scaling law is a 
problem open to further investigation. 

An additional open problem, which we reserve for future studies, is the case of an arbitrary 
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direction of the electric field. It is known that Bloch oscillations of a quantum particle on 
a 2D lattice changes drastically if the field vector is misaligned with respect the lattice axis 



23|. Thus one may expect the quantum transport in the presence of a magnetic field to be 



also different. 
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